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Abstract 

We present a theoretical study of time-dependent quantum transport in a resonant tunnel junc- 
tion coupled to a nanomechanical oscillator within the non-equilibrium Green's function technique. 
An arbitrary voltage is applied to the tunnel junction and electrons in the leads are considered to be 
at zero temperature. The transient and the steady state behavior of the system is considered here 
in order to explore the quantum dynamics of the oscillator as a function of time. The properties 
of the phonon distribution of the nanomechnical oscillator strongly coupled to the electrons on the 
dot are investigated using a non-perturbative approach. We consider both the energy transferred 
from the electrons to the oscillator and the Fano factor as a function of time. We discuss the 
quantum dynamics of the nanomechanical oscillator in terms of pure and mixed states. We have 
found a significant difference between a quantum and a classical oscillator. In particular, the energy 
of a classical oscillator will always be dissipated by the electrons whereas the quantum oscillator 
remains in an excited state. This will provide useful insight for the design of experiments aimed 
at studying the quantum behavior of an oscillator. 
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I. INTRODUCTION 



Nanoscopic physics has been a subject of increasing experimental and theoretical interest 
for its potential applications in nanoelectromechanical systems (NEMS)^^-. The physical 
properties of these devices are of crucial importance in improving our understanding of the 
fundamental science in this area including many-body phenomenal^. One of the most striking 
paradigms exhibiting many body effects in mesoscopic science is quantum transport through 
single electronic levels in quantum dots and single molecules^"- coupled to external leads. 
Realizations of these systems have been obtained using semiconductor beams coupled to sin- 
gle electron transistors (SET's) and superconducting single electron transistors (SSET's)^*^, 
carbon nanotubes^^ and, most recently, suspended graphene sheets^. Such systems can be 
used as a direct measure of small displacements, forces and mass in the quantum regime. 
The quantum transport properties of these systems require extremely sensitive measurement 
that can be achieved by using SET's, or a resonant tunnel junction, and SSET's. In this 
context, NEMS are not only interesting devices studied for ultrasensitive transducers but 
also because they are expected to exhibit several exclusive features of transport phenomena 
such as avalanche-like transport and shuttling instability^'^. The nanomechanical prop- 
erties of a resonant tunnel junction coupled to an oscillator— or a SET^®'^^ coupled to an 
oscillator are currently playing a vital role in enhancing the understanding of NEMS. 

The nanomechanical oscillator coupled to a resonant tunnel junction or SET is a close 
analogue of a molecule being used as a sensor whose sensitivity has reached the quantum 
limiti~-i2i'i^. The signature of quantum states has been predicted for the nanomechanical os- 
cillator coupled to the SET's2 and SSET's^SdE. In these experiments, it has been confirmed 
that the nanomechanical oscillator is strongly affected by the electron transport in the cir- 
cumstances where we are also trying to explore the quantum regime of NEMS. In this system, 
electrons tunnel from one of the leads to the isolated conductor and then to the other lead. 
Phonon assisted tunneling of non-resonant systems has mostly been shown by experiments 
on inelastic tunneling spectroscopy (ITS). With the advancement of modern technology, 
as compared to ITS, scanning tunneling spectroscopy (STS) and scanning tunneling mi- 
croscopy (STM) have proved more valuable tools for the investigation and characterization 
of molecular systems^*^ in the conduction regime. In STS experiments, significant signatures 
of the strong electron-phonon interaction have been observed^"^ beyond the established 
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perturbation theory. Hence, a theory beyond master equation approach or hnear response 
is necessary. Most of the theoretical work on transport in NEMS has been done within the 
scattering theory approach (Landauer) but it disregards the contacts and their effects on the 
scattering channel as well as effect of electrons and phonons on each other^. Very recently, 
the non-equilibrium Green's function (NEGF) approacb^^"— has been growing in importance 
in the quantum transport of nanomechanical systems^^"— i^Ii^^. An advantage of this method 
is that it treats the infinitely extended reservoirs in an exact way^^, which may lead to a 
better understanding of the essential features of NEMS. NEGF has been applied in the study 
of shot noise in chain models^ and disordered junctions^^ while noise in Coulomb blockade 
Josephson junctions has been discussed within a phase correlation theory approach^^. In 
the case of an inelastic resonant tunneling structure, in which strong electron-phonon cou- 
pling is often considered, a very strong source-drain voltage is expected for which coherent 
electron transport in molecular devices has been considered by some workers^^ within the 
scattering theory approach. Inelastic effects on the transport properties have been studied 
in connection with NEMS and substantial work on this issue has been done, again within 
the scattering theory approach^. Recently, phonon assisted resonant tunneling conductance 
has been discussed within the NEGF technique at zero temperature'^^. To the best of our 
knowledge, in all these studies, time-dependent quantum transport properties of a resonant 
tunnel junction coupled to a nanomechanical oscillator have not been discussed so far. The 
development of time-dependent quantum transport for the treatment of nonequilibrium sys- 
tem with phononic as well as Fermionic degree of freedom has remained a challenge since 
the 1980's^. Generally, time- dependent transport properties of mesoscopic systems without 
nanomechanical oscillator have been reported^ and, in particular, sudden joining of the 
leads with quantum dot molecule have been investigated^^'^l for the case of a noninteracting 
quantum dot and for a weakly Coulomb interacting molecular system. Strongly interacting 
systems in the Kondo regime have been investigated^"^. More recently^, the transient ef- 
fects occurring in a molecular quantum dot described by an Anderson- Holstein Hamiltonian 
has been discussed. To this end, we present the following study. 

In the present work, we shall investigate the time evolution of a quantum dot coupled 
to a single vibrational mode as a reaction to a sudden joining to the leads. We employ 
the non-equilibrium Green's function method in order to discuss the transient and steady 
state dynamics of NEMS. This is a fully quantum mechanical formulation whose basic ap- 
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proximations are very transparent, as the technique has already been used to study trans- 
port properties in a wide range of systems. In our calculation inclusion of the oscillator 
is not perturbative as the STS experiments^ii^^ are beyond the perturbation theory. So a 
non-perturbative approach is required beyond the quantum master equation21i2^iii or lin- 
ear response. Hence, our work provides an exact analytical solution to the current-voltage 
characteristics, including coupling of leads with the system, very small chemical potential 
difference and both the right and left Fermi level response regimes. For simplicity, we use 
the wide-band approximation^^'^^i^'^, where the density of states in the leads and hence the 
coupling between the leads and the dot is taken to be independent of energy. Although the 
method we are using does not rely on this approximation. This provides a way to perform 
transient transport calculations from first principles while retaining the essential physics of 
the electronic structure of the dot and the leads. Another advantage of this method is that 
it treats the infinitely extended reservoirs in an exact way in the present system, which 
may give a better understanding of the essential features of NEMS in a more appropriate 
quantum mechanical picture. 



II. MODEL CALCULATIONS 



We consider a single quantum dot connected to two identical metallic leads. A single 
oscillator is coupled to the electrons on the dot and the applied gate voltage is used to tune 
the single level of the dot. In the present system, we neglect the spin degree of freedom and 
electron-electron interaction effects and consider the simplest possible model system. We 
also neglect the effects of finite electron temperature of the lead reservoirs and damping of 
the oscillator. Our model consists of the individual entities such as the single quantum dot 
and the left and right leads in their ground states at zero temperature. The Hamiltonian of 
our simple systent^'^'^ is 

Hdot-ph = [eo + Xl{b^ + b)] cJcq + nw{b^b + i) , (1) 

where eo is the single energy level of electrons on the dot with Cq, Cq the corresponding 
creation and annihilation operators, the coupling strength, rj = XI, with A = eE, is the 
electrostatic field between electrons on the dot and an oscillator, seen by the electrons due 
to the charge on the oscillator, / = a/ h/2mu is the zero point amplitude of the oscillator. 
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u is the frequency of the oscillator and h\h are the raising and lowering operator of the 
phonons. The remaining elements of the Hamiltonian are 



Hxe^As = ^^jc]cj, (2) 
3 

-f^leads-dot = ^ Va{t) {c]cf) + cjcj j , (3) 

j 

where we include time-dependent hopping Va{t) to enable us to connect the leads a = L,R 
to the dot at a finite time. For the time-dependent dynamics, we shall focus on sudden 
joining of the leads to the dot at t = 0, which means Va{t) = V6{t), where 6{t) is the 
Heaviside unit step function. is the total number of states in the lead, and j represents 
the channels in one of the leads. For the second lead the Hamiltonian can be written in the 
same way. The total Hamiltonian of the system is thus H = i^dot-ph + -f^icads + -f^ieads-dot • We 
write the eigenfunctions of -ffdot-ph as 

*„(ir,Xo ^0) = A^exp[-^]H^(/ir)exp[-iirxo] (4) 
^„(ir,xo = 0) = A„exp[-^]H„(/X), (5) 

for the occupied, Xq 7^ and unoccupied, Xq = 0, dot respectively, where Xq = X/2mu'^ 
is the shift of the oscillator due to the coupling to the electrons on the dot, where An = 
1 / a/\7^2^^77IZ, Am = 1 / a/a7^2™7ttT/ , and B.n{iK) are the usual Hermite polynomials. Here 
we have used the fact that the harmonic oscillator eigenfunctions have the same form in 
both real and Fourier space {K). 

In order to transform between the representations for the occupied and unoccupied dot 
we require the matrix with elements = / ^^(-^5^0 = 0)'^rn{K,xo 7^ 0) dK, which may 
be simplified^ as 

<l>n,m = , ! , , / exp {-I'K') H:(/i^) Rm{lK) exp {iKxo) dK (6) 

for n < m, where x = Xq/1 and L^~"(x) are the associated Laguerre polynomials. Note that 
the integrand is symmetric in m and n but the integral is only valid for n < m. Clearly the 
result for > m is obtained by exchanging m and n in equation ([71) to obtain 



l2\m-n\ j^jn 



In order to calculate the analytical solutions and to discuss the numerical results of the 
transient and steady state dynamics of the nanomechanical systems, our focus in this section 
is to derive an analytical relation for the time dependent effective self-energy and the Green's 
functions. In obtaining these results we use the wide-band approximation only for simplicity, 
although the method we are using does not rely on this approximation, where the retarded 
self-energy of the dot due to each lead is given hy^'^ 

Kitut2) = v:ih)g:^jh,h)v.ih), (9) 

where a = L,R represent the left and right leads and the Green's function in the leads for 
the uncoupled system is 

+00 

sL,a{ti,t2) = jr'^gljituh) = -mj{ti -t2) I d£„exp[-i£„(ti -^2)], 
■' —00 
+00 

w,th the fact that - \ ^-^'^^^^ where , stauds for every channel m each lead a,rd „„ 

—00 

is the constant number density of the leads. 

Now using the uncoupled Green's function into equation ([2]), the retarded self energy may 
be written as 



+00 

^^(^1,^2) = -in,^(ti - ^2) y d£«V;*(ti) exp[-i5,(ti - hW^ih), (10) 

—00 

+00 

= -\naV*{tx)Va_{t2)Q{tx-t2) j de„exp[-ie„(ti - ts)], 

—00 

= -mo.V:{h)Vo,ih)e{t, - h) X 27r(5(ti - h), (11) 

Now we use the fact that Va(ti) = \V\ x d{ti), 14(^2) = \V\ x d{t2). Then the above 
expression can be written as 

K{tiM) = -\^TJ{t2)5{h-t2) (12) 
where Ta = AT:\V\^na is the damping factor (T l = Tr = T). Similarly S^()f:i,t2) = 

K(ti,t2)]* = +|ira e{t2)6{h~h) . 
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We solve Dyson's equation using -ffdot-ieads, as a perturbation. In the presence of the 
oscillator, the retarded and advanced Green's functions on the dot, with the phonon states 
in the representation of the unoccupied dot, may be written as 

m k 

where g^"^) is the retarded (advanced) Green's function on the occupied dot coupled to the 
leads may be written as, 

gl,{tM) = -^0{t-ti)^eM-i{e^-iT){t-h)l ti>0 (14) 

g^(t2,t') = +1^(^2-0 X exp[-i(£fc + ir)(t2-t')]> ^2>0 (15) 

with Em = + {m + \)huj — A, = eo + (A; + ]^)f\uj — A, and A = A^/2ma;^. 

The above Eqs. ( JT2l [T3l [T4l [T5l) will be the starting point of our examination of the time- 
dependent response of the coupled system. These functions are the essential ingredients 
for theoretical considerations of such diverse problems as low and high voltage, coupling of 
electron and phonons, transient and steady state phenomena. 

III. TIME-DEPENDENT DOT POPULATION p{t) 

The density matrix is related to the dot population through pit) = ^p„,^(t,t), where 

n 

the density matrix p„,„(t,t) = —iG^^,{t,t'), for t = t' and n = n'. G^„/(t,t') is the lesser 
Green's function2ii2^i2^ on the dot including all the contributions from the leads. The lesser 
Green's function for the dot in the presence of the nanomechanical oscillator is given by 

G<x(^'0= E //dMt2G;,„„(t,ti)S< „, Jti,t2)G:,,„,(t2/), tandt'>0 (16) 

no,nQ,a 

whereas, for t and t' < 0, the G,^„/(t, if:') is equal to zero, and G^^,(t, t') includes all the infor- 
mation of the nanomechanical oscillator and electronic leads of the system, and no,n'f),n,n' 
are the oscillator indices. The lesser self-energy, , (^1,^2), contains electronic and os- 
cillator contributions. The electronic contributions are non-zero only when ti and ^2 > 0. 
As the oscillator is initially in its ground state, only the hq = tiq = term gives a non-zero 
contribution to the lesser self-energy. The lesser self-energy for the dot may be written as 
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with 



+00 

Sa,aitl,t2) = J^^&ajitl^h) = / d^^f « (£a)2m„ exp [-k^. (tl -^2)], 

—00 

where iai^a) is the Fermi distribution functions of the left and right leads, which have 
different chemical potentials under a voltage bias. For the present case of zero temperature 
the lesser self-energy may be recast in terms of the Heaviside step function 6{x) as 

+00 

^ofijtu h) = iFa J {tFa + l^uj - So) ^(^1)^(^2) exp[-i£,(ti - ta)] , (17) 

—00 

where Sq Q"^'*-^''(ti, are all non-zero only when both times (^1,^2) are positive ^1,^2 > 
and epo is the Fermi energy on each of leads. 

The density matrix Pn,n(t, t) can be calculated by using Eqs. ( fT2l [T3l [HI [151 [E]) in Eq. ( !T6|l 
at t = t' and n = n' as 



Pn,n{t^t) = -^y2 f f dtldt2$n,m*S,meXp[-i(£m " i^){t - h)] 
„^^uJO Jo 



a,m,k 



x{ir / ^exp[-i£,(ti-t2)}$o,fc$;,fcexp[-i(£fc + ir)(t2-t)], 



00 



Although g'^'('*)(ti,t2) is non-zero for t < 0, it is never required due to the way it combines 
with Sq Q'^^'*-^''(ti, ^2)- By carrying out the time integrations, the resulting expression is 
written as 

"T , r<^Fa S) <T)* <T)„ , (b* 



X {1 + exp[i(efc - + 2ir)t] - exp[-i(eQ - Ek - iT)t\ - exp[i(£„ - + ir)t]} 

The integral over the energy in the above equation is carried out^^. The final result for the 
density matrix is written as 

ZTT ^ Ek- Em + 2ll 

m,k 

where we have added the contribution from the right and the left leads, which can be written 
in terms of a as 
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F^fc = {1 + exp[i(£fc -em + 2iT)t]} {ln(eFa - £fc - iF) - ln(eFa -e^ + if)} 



{l + exp[i(£fc-e„ + 2ir)t]} 



, ln[(eF^-e,)^ + n 
2ln[(eF,-£^)2 + r2] 



+ i 



r 



r 



+ 7r 



Z";, = exp[i(£fc - + 2ir)t] {- Ei[i(eFa - £fc - ir)t] + Ei[-i(eFa - + ir)t]} 
+ {Ei[i(eFa -em + iT)t\ - Ei[-i(eFa -Sk- iT)t\] , 

with e-poi being the right and the left Fermi levels and Ei(x) the exponential integral function. 
Special care is required in evaluating the Ei(x) to choose the correct Riemann sheets in order 
to make sure that these functions are consistent with the initial conditions p(0) = and 
are continuous functions of time and chemical potential. The same applies to complex 
logarithms in the first, apparently simpler, form for V^^. 

Now using equation (fTSi) . the dot population may be written as 



r 



2tx ^ 



Em + 2ir 



IV. TIME-DEPENDENT CURRENT FROM LEAD a 

The particle current Jq, into the interacting region from the lead is related to the expec- 
tation value of the time derivative of the number operator = '^j'^aj, as^^i^"— 



— le 



{[H,x]) 



(19) 



and the final result for the current through each of the leads is written as (See appendix lAj) 



where 



rlo 



T2a 



Ia{t) 



eV 



2 ( tan"^ 



r 



TT 

+ 2 



-i {Ei[+i(eFa -em + ir)t] - Ei[-i(eFa - £m - ir)t]} 



;i + exp[-2n]) tan 



-1 



-iiexp[-2n] {Ei[+i(eF„ - - ir)t] - Ei[-i(eFa - £™ + iT)t]} 
+ ii {Ei[+i(eFa -em + iT)t] - Ei[-i(eFa - Sm - ir)t]} , 



(20) 



where in calculating the left current we need and both the contributions and and 
for the right current is replaced by I^. As before, special care is required in evaluating 
the Ei(x) to choose the correct Riemann sheets in order to make sure that these functions 
are consistent with the initial conditions Ia{t) = and are continuous functions of time and 
chemical potential. 



V. AVERAGE ENERGY AND FANO FACTOR 



To calculate the energy transferred from the electrons to the nanomechanical oscillator, 
we return to the density matrix Pn,n(t, t) given in Eq. ( |T8i) . We may therefore use the lesser 
Green's function or density matrix to calculate the energy transferred to the oscillator as 

Y,nhujpn,n{t,t) 

E,n = {nhw) = Vn ■ (21) 

n 

Note that the normalisation in equation (12T|) is required as the bare density matrix contians 
both electronic and oscillator contributions. The trace eliminates the oscillator part, leving 
the electronic part. In order to further characterize the state of the nanomechanical oscillator 
we investigate the Fano factor for the change of average occupation number, < n > as a 
function of time. The corresponding relation for the Fano factor is given by^ 

(22) 

where (n) = Yj^Pn,rSt,^MY^Pn,n{t,t) and {r?) = Y.^'^ Pn,n{t,t) /Y.Pn,n{^,t), with the average 

n n n n 

evaluated using the diagonal element of the density matrix on the quantum dot. 



VI. DISCUSSION OF RESULTS 



The dot population, net current through the system, total current into the system, average 
energy and Fano factor of a resonant tunnel junction coupled to a nanomechanical oscillator 
are shown graphically as a function of time for different values of coupling strength, tunneling 
rate, and voltage bias. The following parameters^"— were employed: the single 
energy level of the dot eo = 0.5, and the characteristic frequency of the oscillator ^ = 0.1. 
These parameters will remain fixed for all further discussions and have same dimension as of 
hw. We are interested in small and large values of tunneling from the leads, different values 
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of the coupling strength between the electrons and the nanomechanical oscillator, and of the 
left chemical potential < ep^ < 1. 




Time 

FIG. 1. Time-dependent dot population p{t) against time for different pairs of the right and the 
left Fermi energies (0,0), (0,1), (1,1). The dotted line correspond to empty, dashed line correspond 
to half full and solid line corresponds to almost full state of the dot. Parameters: eo = 0.5, ^ = 
0.1, r = 0.1,77 = 0.05 . Units: all the parameters have same dimension as of Hlj. 




Time 

FIG. 2. Total current {lL{t) + lR{t)) flowing onto the dot as a function of time for fixed values 
of eo = 0.5, hio = 0.1, r = 0.1, r] = 0.05, epR = 0, e^L = l.This current (solid line) is equivalent to 
the rate of change of dot population ^p{t) (dashed line) as a function of time for same parameters 
as of current. In this figure, solid and dashed lines have same values at all points. Units: all the 
parameters have same dimension as of 

The nanomechanical oscillator induced resonance effects are clearly visible in the numer- 
ical results. It must be noted that we have obtained these results in the regime of both 
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strong and zero or weak coupling between the nanomechanical oscillator and the electrons 
on the dot. The tunneling of electrons between the leads and the dot is considered to be 
symmetric (F/j = F^) and we assume that the leads have constant density of states. 

The dot population is shown in fig. [T], as a function of time in order to see the transient 
and steady state dynamics of the system. We consider here empty, half full and occupied 
states of the system for fixed values of F = 0.1, t] = 0.05, by choosing the right and the left 
Fermi levels pairs ( 0, 0), (0, 1) and (1, 1) respectively. Firstly, when both the Fermi levels 
are below the dot energy then the dot population rises initially for a short time and for long 
times settles at a small but finite value. This is not quite empty because the finite F allows 
some tunneling onto the dot. Secondly, when the left Fermi level is above the dot energy 
then the dot population settles in a partially full (half full) state. Thirdly, when both the 
Fermi levels are above the dot energy, it is completely full for a short time but for long time 
is not quite full, again due to the dot coupling with the leads. These results are consistent 
with the particle-hole symmetry of the system as the empty state of the system is not empty 
and the occupied state is not completely full, while the partially full is roughly half full. 

In fig. 121 we have shown the total current fiowing onto the dot as a function of time for fixed 
values of F = 0.1, ?7 = 0.05, epi? = 0, and of the left Fermi level 1. This current (solid line) is 
equivalent to the rate of change of the dot population (dashed line) for the same parameters. 
In this figure, we can not distinguish the solid and the dashed line. This confirms that our 
analytical results are consistent with the equation of continuity, Iiit) + /i?(t) = ^p(t), and 
hence, with the conservation laws for all parameters. 

In fig. |3] we have shown the net current (Iiit) — Init)) fiowing through the system as 
a function of both time and of the left Fermi level for two different values of coupling 
strength: rj = 0.02 to rj = 0.08 and for small and large values of F. We observe simple 
oscillations in the net current flowing through the system for weak coupling strength and 
weak tunneling. With increasing coupling strength the structure of the oscillations becomes 
more complicated as shown in flg. [3t^b). In order to interpret this complicated structure, we 
have a two step discussion: flrstly, we have plotted the net current as a function of time in 
flg. IHwith flxed values of the Fermi level, epL = 1? ^fr = 0, tunneling energy, F = 0.01 and 
for different values of coupling strength: rj = 0.02 and f] = 0.08. In this flgure, in the limit of 
weak coupling the oscillations are again simple while for the strong coupling limit, there is 
a beating pattern in the oscillations. We note that the frequency of the simple oscillations 
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Net Current 



Net Current 




FIG. 3. Net current {I hit) — iRit)) flowing through the system as a function of both time and 
of the left Fermi level for two different values of coupling strength: ij = 0.02 (Fig. EJ^a)), and 0.1 
(Fig. El^b)). Parameters: eo = 0.5, epR = 0,€fl = l,huj = 0.1, F = 0.01. Units: all the parameters 
have same dimension as of hui. 

is (Icfl — Col) and these oscillations are present even in the limit of weak coupling. We 
conclude that this is a purely electronic process (plasmon oscillations). It is clear from the 
figure that in the strong coupling case, it contains two beating frequencies, therefore we 
interpret this as due to a mixture of electronic and mechanical frequencies. Secondly, in 
fig. 13 we have plotted the net current for fixed values of epL = 1, = 0, tunneling energy, 
r = fiw and for different values of coupling strength: rj = 0.02 and rj = 0.08. We have found 
that in the regime ( F > fiui), the effects of the oscillator are not apparent and the period of 
the nanomechanical oscillator can not be resolved. Why can the period of the oscillator not 
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FIG. 4. Net current {I lit) — iRit)) flowing through the system as a function of time for two 
different values of coupling strength: rj = 0.02 (dotted line), and 0.08 (solid line). Parameters: 
eo = 0.5, epR = 0,eFL = l,fiuj = 0.1, F = 0.01. Units: all the parameters have same dimension as 
of huj. 
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FIG. 5. Net current (Iiit) — iRit)) flowing through the system as a function of time for two 
different values of coupling strength: r] = 0.02 (dotted line), and 0.08 (solid line), and F = 0.1. All 
the parameters are same as in fig. |3]and have same dimension as of fuj. 

be resolved by the electrons in this limit? In this regime, electrons spend less time on the 
dot than the period of the oscillator. Therefore, electrons do not resolve the period of the 
nanomechanical oscillator. Now we will focus only in the regime of small tunneling F < hw, 
for further discussion in order to analyze the dynamics of the nanomechanical oscillator and 
the effects of coupling between the electrons and the nanomechanical oscillator. 

Next we have shown the average energy of the nanomechanical oscillator as a function of 





14 



FIG. 6. Average energy transferred to the oscillator as a function of time and left Fermi level for 
fixed values of eo = 0.5, epR = 0, P = 0.01 and for different values of coupling strength: rj = 0.02 
(Fig. Ola)), and 0.08 (Fig. EKb)). Units: all the parameters have same dimension as of fkj. 

time and of the left Fermi energy in fig. [6]for fixed values of tunneling F = 0.01, ep/j = 0, and 
for different values of coupling strength 1] = 0.02, rj = 0.08. We found damped oscillations 
for short times and constant energy for long times. This constant average energy increases 
with increasing Fermi level. Why have we found this particular type of structure? We know 
that the nanomechanical oscillator potential seen by the electrons on the dot is independent 
of time when the oscillator is in any of its pure eigenstates. Otherwise, when the oscillator 
is not in a pure state, the potential seen by the electrons is time dependent. In the former 
case, the electrons are scattered elastically by the time independent potential and in the 
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latter case the scattering process is inelastic because the time dependent potential allow the 
transfer of energy between the two. We observe that the constant average energy also has 
steps as a function of the left Fermi level which become more pronounced with increasing 
coupling strength. Hence, the oscillatory part of the behavior of the mechanical oscillator is 
damped by coupling with the electrons on the dot but the constant part is not. The damping 
mechanism in the transient dynamics is due to transfer of energy from the nanomechanical 
oscillator to the electrons on the dot while when the oscillator is in any of the pure eigenstate 
then there is no mechanism for the transfer of energy between the two. This same physical 
phenomenon also applies to the net current flowing through the dot as well. This appear to 
be a specifically new quantum phenomena in the study of nanomechanical systems. 




FIG. 7. Average energy transferred to the oscillator as a function of j- and for fixed values of 
eo = 0.5, t = WOO,eFR = 0,€fl = l,r = 1 and rj = 0.02. Units: all the parameters have same 
dimension as hw. 

Can we compare this quantum phenomena with the classical mechanical oscillator? Yes, 
the nanomechanical oscillator has to enter the classical regime in the limit of small h. For 
this, we study the dynamics of the quantum oscillator in the classical limit, in which h in the 
mechanical oscillator part of the Hamiltonian given in Eq. ([1]) goes to zero, where hw <T . 
To see this, we have plotted the average energy as a function of y = in the nanomechanical 
part of the system in fig. [7] for fixed values of tunneling F = 1, epij = 0, epL = 1 and coupling 
strength rj = 0.05. We found that the average energy of the quantum nanomechanical 
oscillator scales as h^. We set the average energy in the limit /i' — )■ to see what happen 
to the system for long time. It implies that in this limit, the energy transferred to the 
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nanomechanical oscillator is zero for long time. Hence, we conclude that the long time 
dynamics of the classical mechanical oscillator is always zero. 




20 40 60 80 100 



Time 

FIG. 8. Fano factor as a function of time for two different values of coupling strength: rj = 0.02 
(dotted line), and 0.08 (solid line). Parameters: eo = 0.5, epR = 0,eFL = l,fiuj = 0.1, F = 0.01. 
Units: all the parameters have same dimension as of tjw. 

Finally, in fig. [HI we have shown the Fano factor as a function of time for two different 
values of ?7 = 0.02, t] = 0.08 and for fixed values of F = 0.01, epij = 0, epL = 1. In the limit 
of weak coupling, the nanomechanical oscillator shows thermal like behavior and poissonian 
statistics while in the limit of strong coupling its dynamics is non-thermal which leads to 
super-poissonian statistics. In this figure, the short time behavior is always thermal, but 
this is trivial as the nanomechanical oscillator is initially in its ground state. 

In conclusion, we have found mixed and pure states in our results which confirm the 
quantum dynamics of our model with the following justifications: in a classical mechanical 
oscillator modeU^"— all states give rise to a time dependent potential. Hence, all states of 
the classical mechanical oscillator are damped. Thus, we confirm the new quantum dynamics 
of the nanomechanical oscillator that will be helpful for further experiments beyond the 
classical limit to develop better understanding of NEMS devices. 

VII. SUMMARY 

In this work, we analyzed the time-dependent quantum transport of a resonant tunnel 
junction coupled to a nanomechanical oscillator by using the nonequilibrium Green's function 
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approach without treating the electron phonon couphng as a perturbation. We have derived 
an expression for the full density matrix or the dot population and discuss it in detail for 
different values of the coupling strength and the tunneling rate. We derive an expression 
for the current to see the effects of the coupling of the electrons to the oscillator on the 
dot and the tunneling rate of electrons to resolve the dynamics of the nanomechanical 
oscillator. This confirms that electrons resolve the dynamics of nanomechanical oscillator 
in the regime Te > tqsc while they do not in the opposite case Te > rosc- Furthermore, we 
discuss the average energy transferred to oscillator as a function of time. We also discuss 
the Fano factor as a function of time, which shows thermal behavior and poissonian to non- 
thermal and super-poissonian behavior. We have found new dynamics of the nanomechanical 
oscillator: pure and mixed states, which are never present in a classical oscillator. These 
results suggest further experiments for NEMS to go beyond the classical dynamics. 

Appendix A: 

The particle current Jq, into the interacting region from the lead is related to the expec- 
tation value of the time derivative of the number operator Na = '^j'^aj, a:^^'— 

m = -e{j^x) = ^{[H,x]) (Al) 
Ut) = |{G<Jt,t)V;,o(t) - V*Jt) G<o(t,t)}, (A2) 

where we have the following relations 

G<„(t,t) = J dt' {Gl,{t,t')VoAt')g^At',t) + G<o{t,t')VoAt')gaAt',t)} (A3) 
G<o(^,^) = / dt' {glJt,t')V^,o{t') G<o(t',t) + g< Jt,t')K,o(t') Gl,{t',t)} , (A4) 

where g'a^a'^'^\t,t') refers to the unperturbed states of the leads and given as 

+00 
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with the fact that ^ ^ / iV„„d.., with being the constant numbe. density of the leads 

— oo 

and other uncoupled Green's function in the leads are 

+00 



—00 

+00 



Sa,ait,t') = ^^ia{£a)9a,jit,t') = j d^a f a (ea)i^a exp [-ie„ (t - t')] , 



Now using equations flA3l & IA4p in the equation (IA2p of current through lead a as 



Ut) = I / dt' Tr { (GS,o(t, g<.(t', t) + G<o(t, t')Vo,.(t') g^,„(t', t)) V^,,{t) 



VoUt) (gL, Jt, t')V^,o{t') G<o(t', t) + g<Jt, t')V^,o{t') Gl,{t', t)) } , (A5) 



Using the fact that Eg [^^^''•^''(t', t) = VQ,^{t')ga^a'^^\t',t)Vafi{t), we can simplify the above 
equation as 

m = y dt'Tr {GS,o(t,t')So<o,a(^''^) + GUt^mo,oAt',t) 

- SS,o,Jt,t')G'o<o(i',i) - S<o,Jt,t')Qo(i',^)]} , (A6) 

where Sq Q"^'''^'*(t, t') are non-zero only when both the times {t,t') are positive t,t' > 0. 
Although g^'(")(t,t') is non-zero for t < 0, it is never required due to the way it combines 
with EQQ'*^'*^^''(t, t'). Here we note that we require g''"'(")(t, t') from Eq. ([T4]&[T5]) for positive 
times only (t > 0). The first integral on right hand side of Eq. (1A6P may be solved by using 
Eq. (da [H & as 

Tr fdt'Gl,{t,t')J:<,Jt',t) 
Jo 

V [de^f dt'$o,™$*^exp[-i(£^-ir)(t-t')]exp[-i£„(t'-t)] 

m J Jo 

—00 

E ^o^-*;,... 'h^ | i-expii(..-^^ + irw | 

V <l>o,^$o m {ln(eFa - £m + ir) - Ei[i(eF„ - + ir)t]} , (A7) 



27r 



27r 

ir 



where the final result is obtained using standard integrals^. We note once again that special 
care is required in evaluating the ln(x) and Ei(x) to choose the correct Riemann sheets in 
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order to make sure that these functions are consistent with the initial conditions and are 
continuous functions of time and chemical potential. This statement will also apply to all 
further discussions. 

The second & third integral on right hand side of Eq. (IA6p are written as 



Tr fdt' {G<,{t,t')J:UJt',t) - SS,o,„(t,t') G<o(t',t)} = iF Tr G<o(t, t). 

This integral can be solved in the same way as for the dot population. The final result is 
written as^ 



irTrG<o(t,t) = |^$^<fo,„^$S,™{-(l + exp[-2n]) (t 



tan 



+ 



+ iiexp[-2n] (- Ei[+i(eFa - Sm ~ ir)t] + Ei[-i(eFa - Sm + ir)t]) 

+ ii (Ei[+i(eFa -em + ir)t] - Ei[-i(eFa -Em- ir)t]) | , (A8) 

and the fourth integral on right hand side of equation flA6p can be solved by using Eq. f[T^ 
[la &[I7]) as 

-Tr j\t'T,<,^^{t,t')Gl,{t\t) 
Jo 

= ^J2 j dt'<l>o,r„$S,^exp[-i(e^ + ir)(t' - t)] exp[-ie„(t - t')] 



-ir 

"27 

-ir 

"27 



1 - exp[-i(£a - Sm - i^)t 



^ %,m%,m {ln(eFa -Sm- ir) - Ei[-i(eFa -Sm- ir)t]} 



(A9) 



Using equations (1A71 IA8I fc lA9!) in Eq. (]A6p . the final expression for the current is written 

as 



W = ^ E ^0,^n,m{ll^ + ^"r^ + C), (AlO) 



where components of current are written as 



rla 



T2a 



2 ( tan"^ 



e-pa ~ Sj 

f 



TT 

+ 2 



-i {Ei[+i(eFa -em + iT)t\ - Ei[-i(eFa -e^- ir)t]}], 
:i + exp[-2n])<'tan-^ - 



TT 

+ 2 



-iiexp[-2n] {Ei[+i(eF„ - - ir)t] - Ei[-i(eFa - £™ + iT)t\} 
+ ii {Ei[+i(eFa -em + iT)t] - Ei[-i(eFa -Sm- ir)t]} , 
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where in calculating the left current we need together with both and whereas 
for the right current is replaced by I^. 
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